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walks. The existence, joint continuity and sharp Holder conditions in the set variable of 
the local times of a d-dimensional a-time fractional Brownian motion X = {X(t),t £ 
R + } defined by X(t) = (Xi(t), • • • ,X d (t)), wh ere t > and X\ , • • • , Xd are independent 
copies of Z , are investigated. Our methods rely on the strong local ' nondeterminism of 
fractional Brownian motion. 
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1 Introduction 



In recent years, iterated Brownian motion and related iterated processes have received much 
research interest. Such iterated processes are connected naturally with partial differential 
equations and have interesting probabilistic and statistical features such as self-similarity, 
non-Markovian dependence structure, non-Gaussian distributions; see [2"l fTU } \TT \ [T2" l [T7 1 [T% j 
[28[ [29[ [33j W2\ and references therein for further information. Inspired by these results, we 
consider a new class of iterated processes called a-time fractional Brownian motion (fBm) 
for < a < 2 and < H < 1. These are obtained by taking a fractional Brownian motion 
of index H and replacing the time parameter with a strictly a-stable Levy process Y . More 
precisely, let W = {W(t),t £ M} be a fractional Brownian motion in R with index H, which 
is a centered, real-valued Gaussian process with covariance function 



and W(0) = a.s. Here and in the sequel, | • | denotes the Euclidean norm. Let Y = 
{Y(t),t > 0} be a real- valued strictly a-stable Levy process, < a < 2, starting from 0; see 
Section 3 for its definition and [9[ [36] for further information. We assume that W and Y 
are independent. Then a real-valued a-time fractional Brownian motion Z = {Z(t),t > 0} 
is defined by 



For a = 2 and H = 1/2, this is the iterated Brownian motion of Burdzy [10]. When 
< a < 2 and H = 1/2, Z is called an a-time Brownian motion by Nane [30]. Aurzada 
and Lifshits [3] and Linde and Shi [27] studied the small deviation problem for real-valued 
a-time Brownian motion. Nane |32j studied laws of the iterated logarithm for a version 
of Z. Moreover, when Y is symmetric, for a = 1, 2 and H = 1/2 these processes have 
connections with partial differential operators as described in [2j [31] . 

More generally, it is easy to verify that the process Z has stationary increments and is 
a self-similar process of index H/a. The latter means that, for every constant c > 0, the 
processes {Z(t) : t > 0} and {c~ H / a Z(ct) : t > 0} have the same finite-dimensional distri- 
butions. Gaussian and stable self-similar processes have been studied extensively in recent 
years; see Samorodnitsky and Taqqu [35], Embrechts and Maejima [20] for further informa- 
tion. The a-time fractional Brownian motions form an important class of non-Markovian 
and non-stable self-similar processes, except in the special case when H = 1/2 and Y is a 
stable subordinate!' [In this case, Z is a symmetric stable Levy process]. As will be shown 
in this paper, they have natural connections to partial differential equations and can arise 
as scaling limit of randomly indexed random walks with dependent jumps. Hence they can 
serve as useful stochastic models in many scientific areas including physics, insurance risk 
theory and communication networks. Moreover, because they are non-Markovian and have 
non-stable distributions, new methods are often needed in order to study their properties. 

When a < 2, the sample function of the a-time fractional Brownian motion Z is not 
continuous and its irregularity is closely related to those of W and Y . One of our motivations 
of this paper is to characterize the irregularity of Z in terms of the parameters H and a. 




Z(t) = W(Y(t)), t > 0. 



(1.1) 
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We do this by studying the existence and regularity of the local times of a-time fractional 
Brownian motion X = {X(t),t > 0} with values in M. d defined by 

X(t) = (Xi (t),... ,X d (t)), (t>0), (1.2) 

where Xj = Wj(Yj(t)) (j = 1, ■ ■ ■ , <£). We assume that W\, • • • , W d are independent copies 
of W, Y\,--- , Yd are independent copies of Y, and {Wj} and {Yj} are independent. We 
will call X = {X(t),t > 0} a d-dimensional a-time fractional Brownian motion. It is clear 
that X is also self-similar of index H/a and has stationary increments. 

The rest of this paper is organized as follows. In Section 2, we study the PDE connections 
of a-time fractional Brownian motions, and prove that they can be obtained as scaling limit 
of randomly indexed random walks with dependent jumps. These results provide some 
analytic and physical interpretations for a-time fractional Brownian motions. 

In Sections 3 and 4 we investigate the existence, joint continuity and sharp Holder 
conditions in the set variable of the local times of a d dimensional a-time fractional Brownian 
motion X. In the special case of d = 1, W is Brownian motion and Y is a symmetric a- 
stable Levy process with a > 1, the existence and joint continuity of the local time of Z 
have been proved by Nane [30] . The methods used in this paper differ from those of Nane 
|30j . The latter uses the existence of local times of Brownian motion in K as well as the 
existence of local time of symmetric stable Levy processes, which does not exist whenever 
a < 1 or d > 1. Our Theorem 13.11 implies that, for the a-time Brownian motion X in M. d , 
local times exist in the case d = 1 for a > 1/2; in the case d = 2 for a > 1 and in the 
case d = 3 for a > 3/2. Moreover, Theorem 14.11 shows that these local times have jointly 
continuous versions. 

The methods of Sections 3 and 4 rely on the Fourier analytic argument of Berman [7J [8] 
and a chaining argument in Ehm p3j. In order to derive crucial moment estimates for the 
local times, we make use of the strong local nondeterminism (SLND) of fractional Brownian 
motion proved by Pitt [34] as well as several nontrivial modifications of the arguments in 
Xiao 0HH2]. 

Finally, we provide some technique lemmas as an Appendix in Section 5. 
Throughout the paper, we will use K to denote an unspecified positive finite constant 
which may not necessarily be the same in each occurrence. 

2 PDE connections and scaling limits of randomly indexed 
random walks 

In this section we show that a-time fractional Brownian motions have natural connections 
to partial differential equations. They may also arise as scaling limit of randomly indexed 
random walks with dependent jumps. These results provide some analytic and physical 
interpretations of a-time fractional Brownian motions. These results show that a-time fBm 
can serve as useful stochastic model in various scientific areas. 
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2.1 PDE connections 

The domain of the infinitesimall generator A of a semigroup T(t) defined on a Banach space 
% is the set of all 99 £ % such that the limit 

lim T(t)^)-^ 



i->0 



exists in the strong norm of %. 

Let A = Y2j=i be the Laplacian operator, and let 5(x) be the Dirac-delta function. 

The density of Brownian motion W in M. d is f(t,x) = e - ' 2 ^ 2 / 2 *. Let 



T(t)ip(x)= f(t,x - y)(p(y)dy 



be the semigroup of Brownian motion on L 2 (M rf ). Then the generator of T{t) is A with the 
domain Dom(A) = {<p € L 2 (M d ) : Vy? € L 2 (M ci )}, where V99 is the weak derivative of if. 
See Section 31 in Sato [36] for more details and semigroups on other Banach spaces. Let 
<p be a function in the domain of the Laplacian. Then the function u(t, x) = T(t)ip(x) is a 
solution of the heat equation 

d 1 

—u(t,x) = -Au(t,x), t>0,xeM. d , 

dt y ' 2 v ' " ' (2.1) 

u(0,x) = (p(x), x £ M. d . 

In the case H = ^ and y(£) is stable subordinator of index /3/2, < /3 < 2 with 

is a symmetric stable process of index /3 in R rf . The density 

of W(Y(t)) is given by 

poo poo g— |x| 2 /2s 

q(t,x) = J^ f(s,x)p t (s)ds = (2 7TS )d/2 Pt ^ ds ' 
where pt(s) is the density of Y(t). Then the function 

/•oo 

u(t,x)=E x [iP(W(Y(t)))}= / [T(sWx)]pt(s)ds 

Jo 

is a solution of 

^u(t,x) = -2-^ 2 {-Af/ 2 u{t,x), t>0, xeR d (2.2) 
u(0,x) = ip{x), x€R d , 

where — (— A)^/ 2 is the fractional Laplacian with Fourier transform 

e i ^ x )[-(-Af^{x)]dx = -\kf f e i{k ' x) i;(x)dx, 
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for functions ip in the domain of the fractional Laplacian, see [361 Theorem 31.5 and Example 
32.6]. 

For the case of H = ~ and a = 2, Allouba and Zheng [2] and DeBlassie [IT] showed that, 
for any function ip in the domain of the Laplacian, the function u(t,x) = ~E x [(p(W (Y (t)))] 
solves the Cauchy problem 



0_ 

01 



u(t,x) = + lA 2 u(t,x), t>0, xeR d 



V8nt 

tt(0,x) = tp(x), x G R d . 
In this case u(t,x) = M x [(p(W (Y (t)))] also solves the fractional Cauchy problem 
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T u(t,x) = 2^ 2 Au(t,x); x£R d , t>0 
it(0,x) = (p(x), x G R d . 



Here ^-u(t, •) is the Caputo fractional derivative with respect to t of order |, defined by 
(for fixed x G R d ) 



1 

— = -f= 

9i2 A 



L 



<9ii(s, x) 



ds 



ds 



1 



S)2 



r(i-i) 



<9u(s, x) 
ds 



ds 



. i ' 

J(t- S )2 



(2.3) 



see @]. 

For the case = i, a = 1 and Y is a symmetric Cauchy process, Nane [31] showed 
that u(t,x) = E x [<p(W(Y(t)))} solves 



d 2 , ^ 



A(p(x) 

TTt 



-A 2 u(t,x), t > 0, xG 



it(0, x) = </?(x), 



x G 



(2.4) 



a 2 y 



where ip is a bounded measurable function in the domain of the Laplacian, with q x .q x . 
bounded and Holder continuous for all 1 < i, j < d. 

Nane [31] has also established pde connection for the case H = |, and a = — for 
relatively prime integers k, m, see Theorem 2.5 in [31] , 



For the case 0<H<1, d = l and a = 2, D'ovidio and Orsingher [18] established the 
fact that the density of Z(t) = W{Y{t)) 



q(t,x) 



is a solution of the first order PDE 



e 2s 



2 s 2 



V / 5W^ y/^rf 



ds 



L dq(t,x) 



H d 



(xq(t,x)), t > 0, x G R. 



(2.5) 



dt 2 dx 

Let pt{s) be the density function of the symmetric Cauchy process Y(t) and let 5(x) 
denote the Dirac-delta function. The following theorem answers a question in [18] and we 
prove it and Theorem 12.21 below for the more general case that If is a fractional Brownian 
motion with values in R d . 
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Theorem 2.1. In the case < H < 1 and a = 1 i/je density function 



q(t,x) 



Jo 



f H (s,x)p t {s)ds 



Jo 



e 2 S 



2 
TFT 



(27rs 2H ) d / 2 ir{t 2 + s 2 ) 



ds 



ofW(Y(t)) solves the PDE 

d 2 q(t,x) 2HI m/2] (H) 



8t 2 



nt 



AS(x) - H(2H - l)AG (2H _ 2)jt g(t, 



H 2 A 2 G {4H _ 2)>t q(t, x), x£R d ,t>0, 



(2.6) 



where 



poo 

G^ t q(t, x) = 2 s^ Pt (s)f H (s, x)ds, 7 + 0, 

JO 

and Go,t is i/ie identity operator. 

An operator similar to the operator G 7> ( was introduced in [22]. We refer to their 
Proposition 3.6 and Remark 3.7 for some nice properties of that operator. For the case 
H 7^ 1, it might be a challenging problem to find the right class of functions eft and establish 
the Cauchy problem that is solved by u(t,x) = E x (cj)(W(Y(t)))). This is due to the fact 
that v(t,x) = IL X ((J)(W (t))) is not a semigroup on a Banach space. The general theory of 
semigroups and their generators will not apply in this case 

Proof. Recall that the density function of symmetric Cauchy process Y(t) is 

t 



Pt(s) 



n{t 2 + S 2 ) 



t>0, s G 



Since 



and for t > 



dt 2 



Pt(s) 



-2t(3s 2 -t 2 ) 
(t 2 + s 2 ) 3 



f H (s,x) 



d 2 



ds 



e 2s 



(27TS 2H ) d / 2 



-2t(3s 2 - t 2 ) 



(t 2 + s 



2\3 



ds < 00, 



we apply the Dominated Convergence Theorem to verify the following interchange of the 
second derivative in t: 



d 2 f°° d 2 

0^l(t,x) = 2 f H (s,x)-^p t (s)ds. 



(2.7) 



By using integration by parts to (|2.7p and the facts 



2 



d 2 \ 



ds 2 dt 2 



\Pt(s) = 0; 



0_ 

ds 



{21 



f H (s,x) = Hs 2H - 1 Af H (s,x) 



we derive 

d 2 



dt 2 



q(t,x) 



d 2 

f H (s,x)-^p t (s)ds 







B oo r°° f) 2 

2 Pt (s)-^-f H (s,x) q -2j Pt ( s )—f H (s,x)ds 

f) OO f'°° f) { \ 

2p t (s)-^-f H (s,x) o + 2 J Pt {s)^-(Hs 2H - l M H {s,x)y s 



2HI (0,l/2] (H) 



TTt 



A5(x) 



-2 J™ p t (s) (h{2H - l)s 2H - 2 Af H (s, x) + H 2 s* H - 2 A 2 f H (s, x)\ ds. 

2HI (0,l/2] (H) 



TTt 

r oo 



-AS{x) 



f'OG POO 

-A2 p t (s)H(2H -l)s 2H - 2 f H (s,x)ds + A 2 2 Pt (s)H 2 s AH - 2 f H \s,x)ds. 
Jo Jo 

the last line follows by the dominated convergence theorem. In the above we have used that 



and that 



limf H (s,x)^-p t ( S ) = 0, 

s— >0 OS 

lim f H (s,x)^-p t (s) = 0, 

s->oo OS 

d 

lim Pt(s)—f H (s,x) = 0, 

s->oo OS 



Ihn p t (s)-^-f H (s,x) = hm -s^A^^x) 

s->0 OS s->0 7TI 



(2.9) 







if H > i 



2 

^A(J(a;) if < H < ±. 



This finishes the proof of (|2.6j) . 



□ 



Remark 2.1. After we submitted our paper, we learned that Beghin et al. |6j/ has established 
that the density W(Y(t)) in the case d = 1, < H < l,a = 1 



q{t,x) = 2 



e 2 S ^ 



V^ 2 ^ K(t 2 + S 2 ) 



ds 



solves 



d 2 1 



8 B 2 
H(H-l)^x-H 2 ^- 2 x 2 
ox ox* 



q(tx) _2jn^majs^ (210) 



Tit 



It is interesting to compare the equations (|2.6|) and (|2.10p . 
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Let < (3 = < 2 for k, m relatively prime integers, and let Y(t) be a stable subordi- 
nates of index (3/2. In this case the density pt(s) is a solution of 

Q2m Qk 

Q^Pt(s) = ^:Pt(s), s,t>0, (2.11) 

see Lemma 3.1 in [16]. We have the next theorem which gives an extension of the PDE in 
(1221). 



Theorem 2.2. Let < H < 1, fj = — for m = 2, 3, • • • . Let Y(t) be a stable subordinator 
of index a = ^. Then the density q(t,x) = f H (s,x)pt(s)ds ofW(Y(t)) is a solution of 

~Q£2m x ) = -HAV [2H _ 1):t q(t, x), (2.12) 

where Vyj q(t, x) = s" f pt(s)f H (s,x)pt(s)ds for 7^0 and Vo,t is the identity operator. 

Proof. The proof follows by integration by parts as in the proof of Theorem |2.1[ and by 
using (12. HI) with k = 1. 

Q2m poo Q2m 



dt 2r 



poo Qlra 



d 



f H (s,x)—p t {s)ds 

OS 

OO POO Q 

= f H (s,x) Pt (s) - / —f H (s,x)p t (s)ds 
o J ds 

= f H (s,x)p t (s)\™ - J™ pt&fHs^AfX&x)) ds 

poo 

= -H / p t (s)s 2H ~ 1 Af H (s,x)ds, 
Jo 

poo 

= -HA / Pt (s)s 2H - 1 f H (s,x)ds, 
Jo 

the last line follows by dominated convergence theorem. See equations (2.7)-(2.10) in [16] 
to show that the boundary terms are all zero. □ 

Letting H = i in Theorem 12.21 the density q(t,x) of symmetric stable process W(Y(t)) 
of index a = — is a solution of 



in 



Q2m 1 

— q(t,x)=--Aq(t,x). (2.13) 

The equations (j2.2[) and f)2. 13|) should be compared. This result is a special case of the 
following result stated in Nane |3H Lemma 3.2] that is due to DeBlassie [16] originally: Let 
H = \ and < a = — < 2, where k, m are relatively prime. Let Y be a stable subordinator 
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of index a/2. In this case W(Y{t)) is a symmetric stable process of index a = —. Then 
the density q(t,x) of W(Y(t)) is a solution of 

Qlm y 

q(t,x) = - ] -(-A) k q{t,x), s,t>0. (2.14) 



dt 2m K ' 2 k 

We can work out a similar connection for the case < H < 1 and a = ^- (k and m are 
relatively prime integers) by using integration by parts, (12.8|) . and (I2.14p . which extends 
the PDE connection in Nane |31| Theorem] for the case H = | , a = . 

Alternatively, for the case < -ff < 1, a = — , k,m relatively prime integers, we have 
W(Y(t)) = W(B(U (t))), where B is a Brownian motion running twice the speed of standard 
Brownian motion and U is a stable subordinator of index a/2 = In this case using the 
methods of Theorem 12. 1\ equation ()2.5[) for the density of W(B(t)) and equation (|2.1ip for 
the density of U(t) we can obtain the PDE solved by the density of W(Y(t)). 

For many other PDE connections of different types of subordinate processes, we refer 

to [2 12 02 [ill EJ [33]. 

2.2 Scaling limits of randomly indexed random walks 

Now we prove that the a-time fBm Z = {W(Y(t)),t > 0}, where W is fl-fractional 
Brownian motion with values in R and {Y(t), t > 0} is a stable subordinator of index 
a G (0, 1), can be approximated weakly in the Skorohod space D ([0, oo), R) by normalized 
partial sums of randomly indexed random walks with dependent jumps. 

Let {£ n , n 6 Z} be a sequence of i.i.d. random variables with E[£] = 0, E[£ 2 ] = 1, and 
let {a n , n G Z + } be a sequence of real numbers such that 

oo 

zZ a l < °°- 

n=0 

We consider the linear stationary process {X n , n G N} defined by 

oo 

X n = Y, a £n-j, nGN. (2.15) 

3=0 

Davydov |15j was the first to study the weak convergence of normalized partial sums of 
{X n , n G N} to fractional Brownian motion. The following result is taken from Whitt |40t 
Theorem 4.6.1]. 

Lemma 2.1. Let {X n , n G N} be the linear stationary process defined by \2.15\) . and let 

S n =X 1 + ---+X n . If 

Var (5„) = n 2H L(n), n G N (2.16) 
for some H G (0, 1), where L(-) is a slowly varying function, and 

E[|S n | 2a ] <K.{E[S 2 n ]) a (2.17) 



S 



for some constants a > 1/H and K > 0, then 

f 1 



n H yfZ[r< 



S [ntj ,t>0\^{W(t),t>0} (2.18) 



in the Ji-topology on D ([0, oo), ]R) , where W is a fractional Brownian motion with Hurst 
index H. 

Example 2.2. As in [40\ pp. 123-124] we take aj = cj~^ for some constants c G K\{0} and 
7 G (^, 1), then it can be verified that 

Var(5 n ) ~ cin 3 " 27 as n -> oo, (2.19) 

where 

_ 2c 2 r(l- 7 )r(2 7 - 1) 

Cl " r( 7 )(3-2 7 )2 ■ 

By applying Lemma 12. 11 we have that 

| -^jS Vntl , t > 1 => {W(t), t>0}, (2.20) 
in the Ji-topology on D ([0, oo), R) , where H = 3 ~» 7 . 

Theorem 2.3. Lei {X n ,n > 1} be the linear stationary process defined by 12.15\) and 
satisfies $2.16\) and \2.11ty in Lemma \2.1\ Let {J n ,n > 1} be a sequence of i.i.d. random 
variables also independent of {£ n , n G Z+}, which belongs to the domain of attraction of 
some stable law Y with index a G (0, 1) and Y > a.s. Denote by {b n ,n > 1} a sequence 
of positive numbers such that b n T n Y , where 

T n = J\ H h J n , Vn > 1. 

Then as c —> oo 

=5 L T(rf)j, t > I t > 0} (2.21) 



(b(c)y H ^L(b(c)-i) 
in the J\-topology on D ([0, oo), R) , where 6(c) = ^cj and T(s) = 

Proof. For i > 0, let 6(t) = 6m. Then 6(i) = i _1 / Q ^(i), for some slowly varying function 
£(t) at infinity. It follows from Theorem 4.5.3 in Whitt [40] that 

{b(c)T(ct), t > 0} =>- {Y(t), t > 0} (2.22) 

in the Ji-topology on Z) ([0, oo), K+) , where {Y(i), i > 0} is a stable subordinator with 
index a. 
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Notice that {X n , n > 1} and { J„} are independent, we derive from Lemma 12.11 and 
(f2722|) that as c -> oo 

|( ^J__ S LetJ , 6(c)T(ci)), i > o| {(W(t), Y(t)), t > 0} (2.23) 

in the Ji-topology on D ([0, 00), R) x D ([0, 00), R+) . 

Since, for the limiting processes in (12.23[) . the sample function W(t) is continuous and 
Y(t) is strictly increasing, the conclusion of Theorem 12.31 follows from Theorem 13.2.2 in 

Whitt na. □ 



3 Existence of local times 

Let X = {X(t),t > 0} be an a-time fractional Brownian motion in R d defined by (|1.2p . In 
this section, we study the existence of local times 

L = {L{x, B): x G R d , B G B(R + )} 

of X, where i3(R+) is the Borel cx-algebra of R + . In Section U we will establish joint 
continuity and sharp Holder conditions in the set variable for the local times. 

We recall briefly the definition of local times. For an extensive survey, see Geman and 
Horowitz [21]. Let X : R -> R d be any Borel function and let B C R be a Borel set. The 
occupation measure of X(t) on B is defined by 

Hb(A) = \x{t G B : X(t) G A} (3.1) 

for all Borel sets A C R d , where Ai is the Lebesgue measure on R. If hb is absolutely 
continuous with respect to the Lebesgue measure on R d , we say that X has a local 
time on B and define its local time L(x,B) to be the Radon-Nikodym derivative of \xb- If 
B = [0, t], we will simply write L(x,B) as L(x,t). If / = [0,T] and L(x,t) is continuous as 
a function of (x, t) G R d x /, then we say that X has a jointly continuous local time on /. 
In this latter case, the set function L(x, •) can be extended to be a finite Borel measure on 
the level set 

Xj l {x) ={tel: X(t) = x}. 

See Adler [lj Theorem 8.6.1]. This fact has been used by many authors to study fractal 
properties of level sets, inverse image and multiple times of stochastic processes. Related 
to our paper, we mention that Xiao [32] and Hu |23] have studied the Hausdorff dimension, 
and exact Hausdorff and packing measure of the level sets of iterated Brownian motion, 
respectively. 

An a-stable Levy process Y = {Y(t),t > 0} with values in R is a stochastically continu- 
ous process with stationary independent increments, Y(0) = 0, and characteristic exponent 
ip given by 

(ITOt \ 
1 — ij3 sgn(z) tan —J , a / 1; 

/ 2 x (3.2) 

ip(z) = —cr\z\ (1 + i— (3 sgn(z) ln(|z|) j , a = 1, 
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where < a < 2, er > and — 1 < j3 < 1 are fixed constants (we have tacitly assumed that 
there is no drift term). See Bertoin [9] and Sato [36] for a systematic accounts on Levy 
processes and stable laws, respectively. 

Throughout this paper, we assume that Y = {Y(t),t > 0} is strictly stable. That is, we 
assume f3 = in (|3.2|) when a = 1 so the asymmetric Cauchy process is excluded. Strictly 
stable Levy processes of index a are (l/a)-self-similar. Recall that, for t > 0, pt{x) is the 
density function of the random variable Y(t). It is a bounded continuous function with the 
following scaling property: 

Pt{x) = p r t{r l I a x) r l I a for every r > 0. (3-3) 

As discovered in Taylor [39], it is natural to distinguish between two types of strictly 
stable processes: those of Type A, and those of Type B. A strictly stable process, Y, is of 
Type A, if 

p t (x) > 0, Vt > 0, x G R; 

all other stable processes are of Type B. Taylor [39] has shown that if a G (0, 1) and Y is of 
Type B, then either Y or — Y is a subordinator, while all other strictly stable processes of 
index a / 1 are of Type A. Hence, without loss of generality, we will assume Y is either a 
strictly stable process of type A, or a subordinator. It will be shown that the properties of 
local times of a-time fractional Brownian motion X in M. d depends on the type of Y. 

The following existence theorem for the local time of X is easily proved by using the 
Fourier analysis (see, e.g., Berman [7], Geman and Horowitz [21] or Kahane [24]). 

Theorem 3.1. Let X = {X(t), t > 0} be an a-time fractional Brownian motion in M. d . 
Then for any T > 0, X has a local time L(x,T) G L 2 (P x A^) almost surely if and only if 
d < a/H. 

Proof. Let /ti^T] De the occupation measure of X on [0, T] defined by (|3.ip . Then its Fourier 
transform can be written as 

fi[o,T\(u) = / exp(i(u,X(t)))dt, 
Jo 

where (■, •) is the ordinary scalar product in M. d . It follows from Fubini's theorem that 

E f \fi{o,T]{u)\ 2 du= [ [ [ Eexp(i(u,X(t) - X(s))) dudsdt. (3.4) 
Ju. d ' Jo Jo Ju d 

To evaluate the characteristic function in (|3,4p . we assume < s < t (the other case is 
similar) and note that the density of (Y(t),Y(s)) is given by 

Ps,t{x, y) = Ps(x)p t -s{y ~ x). 
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Since X±, • • • , are independent copies of Z = {W(Y(t)),t > 0}, we have 

d 



(3.5) 



Eexp(i(u,X(t) - X(s))) = J] Eexp (iu k (W(Y(t)) - W(Y(s)))) 

k=l 

= TT / / Eexp (iu k (W{y)-W(x)))p s , t {x,y)dxdy 

£Ji JR JR 
d f f 2 

= I1 / / exp ( ~ ~i~\ y ~ X \ 2H ) PsA x >y) dxd y 

=n / ^(-I^M^- 

To evaluate the integrals with respect to u, we make a change of variables to get 

/ exp ( - ^ |z| 2H ) du fc = / exp ( - ^) du fc . (3.6) 

JR V ^ 7 JM V ^ 7 



It follows from ([331), (GLSJ) and Q3j>J) that 



In the above, we have used the fact that P|t_ s i(£) = \t — s\~ l l a p\(\t — s\~ 1 / a z) and another 
change of variables. 

The last integral in (|3,7p is finite if and only if dH/a < 1. Hence //(•) £ L 2 (P x A^) if 
and only if dH/a < 1. Therefore, Theorem 13.11 follows from Plancherel's theorem (see also 
Theorem 21.9 in Geman and Horowitz [21]). □ 



We can express the local time L(t, x) as the inverse Fourier transform of jU[o,T] ( u )i namely 

v 7;/;; , < 3 - 8 > 

exp ( — i(u, x)) exp (i(w, X(s))) du ds. 



, 27r / Jo JR d 

It follows from (|3.8[) that for any x, w; £ M^, -B 6 i3(]R + ) and all integers n > 1, we have 
E[L(x,S)]« = (2vr)-" rf / / exp ( - * f>„ x>) 

JB n JR" d \ ■-, / 

(3.9) 



Eexp U^2(uj,X(tj)) I dudt 
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and for all even integers n > 2, 

E[L(x + w,B) - L(x,B)] n 



(2tt) nd / T ( exp(— i(uj,x + w)) - exp(-i(uj,x}) 



(3.10) 



E exp 

(uj,X(tj)) ) dudt, 

i=i 



where u = (in, • • • , u n ), i = (ii, • ■ • ,i n ) and each Uj G M. d , tj E B (j = 1, ■ ■ ■ , n). In the 
coordinate notation we then write Uj = (uj, • • • ,u d ). For details in deriving the equations 
(|3.9p and (|3.10p . see Geman and Horowitz |21| . 



4 Joint continuity and Holder conditions 

In this section, we establish the joint continuity and sharp Holder conditions in the set 
variable for the local times of d-dimensional a-time fractional Brownian motion X. Then 
we apply these results to study the irregularities of the sample paths of X(t). 

We use methods which are similar to those in Ehm [19] and Xiao [41[ [4*2] . The following 
Lemmas 14.11 14.21 and 14.31 give the crucial estimates for the moments of the local time of 
a-time fractional Brownian motion. Note that the estimates in the case Y is of type A [i.e., 
(|4.2p and (14. 3p ] are different from the case when Y is a stable subordinator [see Lemma [4.3j . 

We need the fact that fractional Brownian motion W satisfies the property of strong 
local nondeterminism (SLND), which was proved by Pitt [34]. More precisely, for any 
i/i,-- - ,y n e R, 

Vax(W(y n )\W( yi ), • • • , W(y n -i)) > K min \y n - Vj \ 2H , uu 

0<j<n—l y ' 

where yo = and K > is an absolute constant. 

Lemma 4.1. Let X = {X(t), t > 0} be a d- dimensional a-time fractional Brownian motion 
with d < a/H for which Y(t) is of type A. For any h > 0, B = [0, h], x G M. d , any integer 
n > 1, we have 

E [L(x, B)] n < K n ( n! )^(i+i/«) j (4.2) 

where K > is a finite constant depending on d, H and a only. 

Proof. Thanks to the strong local nondeterminism (SLND) of fractional Brownian motion 
[cf. Eq. (|4,ip j. Lemma l5?2l and the scaling property of pt(x) [cf. Eq. (|3,3p ]. the proof of 
Lemma [4J] follows along a similar line of the proof of Eq. (2.11) in Xiao [42] with obvious 
modifications. We omit the details. □ 

Lemma 4.2. Under the conditions of Lemma \4.1\ we have that for all even integers n>2 
and < 7 < \ mva{a/{Hd) - 1, 1 - H}, x,w £R d 

E[L{x + w,B) -L(x,B)] n < K n \w\ n ^h n{l -^ )H/a \n\) d+ ~i +m ^ , (4.3) 
where K > is a finite constant depending on d, H , 7 and a only. 
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Proof. Even though the arguments for proving Lemma [4. 2 1 are similar to that of Eq. (2.12) 
in Xiao |42j . several essential modifications are needed. 
By (|3.10p and the elementary inequality 

| e *« _ i| < 2 1 ~ 7 |u| 7 for all u G R and < 7 < 1, 

we see that for any even integer n > 2 and any < 7 < 1, 

E[L(x + w,B) - L{x,B)] n 

(4.4) 



: \wH 



I I f\\ u 3\^^(i^y {uj,X(tj))\ dudt. 



By making the change of variables tj = hsj, j = 1, • • • , n and Uj = h~ H / a Vj, j = 1, • • • , n 
and changing the letters s, v back to t, u, the self-similarity of X implies that the right-hand 
side of (|4.4p equals 

^nj h n(l-(d+j)H/a) f f f\ \ U \1 E exp ( if" ( Uj , X (tj))) dudt. (4 5) 

We fix any distinct t\, ■ ■ ■ ,t n G [0, 1] satisfying 

= t < h < t 2 < ■ ■ ■ < t n , (4.6) 

and consider the inside integral in (|4.5p . Since for any < 7 < 1, \a + 6| 7 < |a| 7 + |6| 7 , we 
have 

n , n 

where the summation ^ is taken over all (A7, • • • , k n ) G {1, • • • , d} n . 

Let us fix a sequence (ki, ■ ■ ■ , k n ) G {1, ■ ■ ■ , (i} n , and consider the integral 

„ n , n \ 

J= TT |4 J | 7 IEexp ( i y"( Uj ,X(tj)) ) du 

„ n / d n \ 

= / |^| 7 Eexp li^^viWifXtitj)) du, 



(4i 



since X(t,) = (Wi(Yi(ij)), ■ ■ ■ ,W d (Y d (tj))). Now, we condition on Y t (tj) = y £j , i 
1, • • • ,d,j = l,--- ,n. By independence of the processes Yi we have that the density of 

(Y e (tj) =y ej : i = 1,-- - ,d,j = 1, • • • ,n) 

is given by 

Ptl,-,tn(2/H>' " ' ,2/ln,2/21, • ' ' ,2/2n,-" , , ' ' " , 2/dra) 
c( n 

=nn**j-*i-ifo#~f'C7-i))- 
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Let t = (ti, ■ ■ ■ , t n ) and y = (y u , ■■■ , yi n ,y 2 i, ■■■ , V2n, ■■■ , Vdl, ■ ■ ■ , Vdn)- By conditioning 
we have 



J 



n I d n \ 

jR" d jR" d - =l y £=1 j=l J 



For any fixed y € R , let 

n I d n \ 

J = / TT \ u f \m exp t V V ulWe( yij )) du 

J**fA \tiU ) 

= I f[ |uf r exp ( -ivar ^E^W ) d «- 

• /R " d 3=1 \ V=l 3=1 / / 

Then by a generalized Holder's inequality and Lemma 15.41 we have 

l/n 



3=1 



lit- I"' 7 



exp ^-^ Var fEE^^wj j du 



nd— 1 yi 

< ( 27r ) 2 / M^e-^dwTT — 

" (det Cov(W e (y tj ), 1 < £ < d, 1 < j < n))V2 7 R 1 1 1 = 1 ^ 



n 1 

(det Cov(Wi(Wi). 1 < * < d, 1 < J < n))V2 II 



if n (n!)' 



n 



n?,i(<letCov(H' f (i, ( j), 1 < j < n))i/2 ^ <7^ 
where 

a 2 kjJ = V a v(W kj (y k . d )\ W e (y e J : t + k or I = kj, i± j) 
= Yai(W kj {y kj:j )\ W kj (y kj:i ) : i = or i^j). 

For any i £ {1, • • • , d} and any y^i, • • • , y^ n , there exists a permutation tt£ of {1, • • • , n} 
such that 

ye,m(i) — y£,7r £ (2) < < Ve,w e (n)- 
Hence, if we write fc,- = ^, then by SLND of fractional Brownian motion (|4.ip . 

ajj = V&T(Wi(y itj )\ W^i) : i = or i/j) 

> K min{|y Wi) - y M -!)\ 2H , \ve,* t (j+i) ~ y^ t {j)\ 2H }- (4.10) 

Hence 

n d n 

n^£*"nn 



3=1 k j>3 1=1 j 



A -\ min{|y £)7r£(j ) - |y^(j+i) ~ ^(i)!}^''"^ 
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where rjej = 1 if kj = I and rjej = otherwise. Note that 



d n 

2 2 ^ = n - 

<=1 J=l 



(4.11) 



Since 



d n 

nn 

1=1 j 

d n 



1 1 

+ 



d n 



1 



^ n n V|^ (J) - y^.Dl^W ' 



where the summation ^ is taken over 2 terms and 5gj 6 {0, 1, 2} and, thanks to (|4.1ip . 

d n d n 



we obtain 



£=i j=i j=i 



n „ d n , 



n 



Now we go back to estimating J. Let 



= {{ye,!, ■ ■ ■ ,ye, n ) ■ ye,n e (i) < ye^p) <••■ < ye^n)}- 



Then 



^<wee n 

{n e } 1=1 



n 



dy 



Fix £ G {1, • • • , d} and a term in , we proceed to estimate the integral 

Ptj-tj-AVij ~ Vi,j-i) 



where dye = dyex ■ ■ ■ dy^^ n . Note that we can write (|4.15p as 



dye, 



f yr Pt* t u)-t* lU) -i(yi,m(j) - yt,* t (j)-i) 
A-nA^ f} x (y eMj) - ye,n e u-i)) H(1+Se ' 3l) 



dye- 



(4.12) 



(4.13) 



(4.14) 



(4.15) 



(4.16) 
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It will be helpful to notice the difference of ye^jfi-i in the numerator and Vt^Aj-i) m the 
denominator. 

Since p t (x) = *~ 1/a Pi(s/t 1/a ), for all t > 0, x € R. Now flgjg) can be written as 

„ / n J ye,* t U)-ye,"tU)-i \ x 

11 ft —* - .w«^ — H"* ( 41? ) 



We now integrate in the order dy £j7re (n), ^ )7r/ ( n _i), ' ■ ' >^,^(1)- 
A change of variables 

ye,n e (n) ~ ye,n e (n)-l = (^(n) ~~ ^(n)-l) 1//Q 2 n 

gives 

' 2/1. 



(^(n) - ^,(n)-l) 1/Q (l/Wn) " ^(n-l)) H(1+ ^ 7) 



Vl,Tv t {n-i) 

pi(z 



(4.18) 



^7rf(n) ^TT£(n)-l) 

where 

_ ye,n e {n-l) — yi,ir e (n)-l 
(*7^(n) — ^(n)-l) 1 ^ 

Since we have assumed < 7 < i(l — if), so #(1 + £^7) < 1 for all £, j. Thus 

— pi(z "L a < 

where iT is a constant independent of s n . This can be verified directly by splitting the 
interval [s n , 00) into [s n ,s n + 1] and [s n + l,oo). 
Continuing this procedure we derive 



^ — Wa+snT) dyt 



fJi (yeMi) ~ Ve,-Kz(j-i)) H{X+h ' jl) 
1 



1 3 

n 

<K n \ 



(4.19) 



Combining this inequality with equation (|4.14|) gives 

d n _. 

J<K«(n\y £ nil ^T^) 

<jf»(nir £ II 1^ ; ; } ( 4 - 2 °) 

Tl,-,-7r d j=l f^- — tj-l) ' I U > 

n 

< K n (n\y +d TT , 



17 



where < e-,- < 2d and Y^j=i e j ^ 2n, thanks to (|4.12p . 
Hence we have shown 

E[L(x + w,B) - L(x,B)] n < K n \w\ n ~ i h n{1 -( d+ ^ H l a \n\) d+ ~ t+l 

di 



(4 21] 

< iiT n |u;| n7 /i n ( 1_ ^ +7 ^/ a )(n!) d+7+1 ' 

n-=ir(i-f(d + 7 ^)) 



r(i + n-E-=if(d + 7^))' 

In the above we use Lemma 15.31 with the fact that 77 + 7£j) < 1 because 7 satisfies 
< 7 < i(a/(#d)-l). It is now clear that (|4.3p follows from (I4.2ip and Stirling's formula. 
This completes the proof. □ 

We have similar bounds for the moments of the local time in the case when Y is a stable 
subordinator. It should be noted that the power of n\ in (14,22j) is different from that in 
Lemma 14.11 This will lead to different forms of laws of the iterated logarithm for the local 
times in the two cases. 

Lemma 4.3. Let X = {X(t), t > 0} be a d- dimensional a-time fractional Brownian motion 
with d < a/H for which Y(t) is a stable subordiantor with index a < 1. For any h > 0, 
B = [0, h], x, w £ M. d , any even integer n > 2 and any < 7 < 9 min{a/(Hd) — 1, 1 — H}, 
we have 

E[L(x,B)] n < K n h^- dH ^ n {nl) dH / a , (4.22) 
E[L(x + w,B) - L(x,B)} n < K n \w\ n ^h n ^ d+ ^ H l a \n\y +m ^ Jl , (4.23) 
where K > is a finite constant depending on d, H , 7 and a only. 

Proof. The proof of (|4.22p is similar to the proof of Lemma l4.1i However, since the sample 
function Y(t) is increasing, for t\,--- ,t n that satisfy (I4.6p . the corresponding yj = Y(tj) 
(j = 1, • • • , n) satisfy yi < ■ ■ ■ < y n , which leads to some clear modifications to the proof. 
Equation (I4.23P follows similarly from the proof of Lemma 14.21 □ 

Now we are ready to prove the joint continuity result of local times. 

Theorem 4.1. If d < a/H , then almost surely X = {X(t), t > 0} has a jointly continuous 
local time L(x, t) (x G R d , t > 0). 

Proof. The proof follows from Lemmas 14.11 14.21 14.31 and Kolmogorov's continuity theorem. 

□ 

Remark 4.1. For d = 1, H = 1/2 and a > 1 this result was proved by Nane 130^ . Theorem 
\4-l\ implies that for d = 1, H = 1/2 and a > 1/2, almost surely X(t) (t > 0) has a jointly 
continuous local time L(x, t) (x G M. d , t > 0). Hence Theorem \4-l\ is an improvement of 
the results in 130\/ and an extension of results in obtained for multidimensional iterated 
Brownian motion. 
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The following tail probability estimates are used in deriving the sharp Holder conditions 
in the set variable of the local times of a-time fractional Brownian motion. 

Lemma 4.4. Suppose Y is not a subordinator. For any A > 0, there exists a finite constant 
A > 0, depending on A, d, H and a only, such that for all r > 0, h > 0, B = [r, r + h], 
x,w eM. d , all < 7 < \ mm{a/(Hd) - 1, 1 - H}, and all u > 

f{l(x + X(t),B) > Ah l - dH l a u dH{l+l/a) } < exp(-Au), (4.24) 

f{\L(x + w + X(t),B)-L(x + X(t),B)\ > A\w\^h 1 -^^ H f a u ci - H ^ 

< exp(-Au), (4.25) 

where C(H, a) = d + j+ H ( d + 2 ^ . 

Proof. Since X = {X(t),t > 0} has stationary increments, i.e., for any r > 0, the processes 
{X(t + r) — X(r),t > 0} and X have the same finite dimensional distributions. Hence 
Lemmas 14.11 and 14.21 can be reformulated as follows: For any r > 0, h > 0, B = [r, r + h], 
x, w € M. d , any even integer n > 2 and any < 7 < \ vam{a / {H d) — 1, 1 — ff}, we have 

E[L(x + X(r),£)] n < r/i^-^'^n!)^ 1 ^/"), (4.26) 
E[L(x + w + X(t),B) - L(x + X( T ),B)} n < K^p/^Md+^A*) 

x(n!) d+7+ ^r^, (4.27) 

where -fC > is a finite constant depending on d, H, 7 and a only. 

Now, Lemma [4.4l is a direct consequence of (|4.26p . (|4.27j) and the Chebyshev's inequality. 

□ 

Lemma 4.5. Suppose Y is a stable subordinator of index a < 1. For any A > 0, there 
exists a finite constant A > 0, depending on X, d, H and a only, such that for all r > 0, 
h > 0, B = [t, t + h], x, w G R d , all < 7 < \ mm{a/(Hd) - 1, 1 - and u > 

F^L{x + X{T),B)>Ah l - dH l a u dH/a ^ < exp(-Au), (4.28) 

p{|L(x + «; + X(t),5) -L(a; + X(r),5)| > A^P^-^^n 15 ^ } 

< exp(-Au), (4.29) 

where D(H, a) = 7 + H(d + 27) . 

The proof of Lemma 14.51 follows the same idea as that in the proof of Lemma 14.4^ with 
an application of Lemma 14.31 We omit it here. 

The next lemma shows that process the real- valued process Z(t) = W(Y(t)) has heavy 
tails as in the case of Y. This might make this process more desirable, since it has heavy tails 
without independence of increments and with the stationarity of the increments. We need 
the following lemma to prove Lemma l4.7l for a two sided estimate of p| sup 0<t<1 \Z(t) \ > it| , 
which will be useful in proving Theorem 14.21 
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Lemma 4.6. Let d= 1, < H < I and < a < 2, and let < a < b then 

Pl\Z(b) -Z(a)\ > u\ 
lim — ^ m j - = C(b - a) 

«-+oc U~ a l H 

for some finite constant C > 0. 

Proof. By using the stationarity of the increments and the self-similarity of W and Y we 
get 

P||Z(6) - Z(a)\ >u} = p{|VF(y(6 - o))| > uj 

= P{(6 - a) H / a \Y{l)\ H \W{\)\ > 

= ¥{(b-a) H / a \Y(l)\ H \s\>u}f H (s)ds (430) 
= [°° f!\Y(1)\ >u llH (b-a)- 1 l a \s\- l l H \f H {s)ds, 

J — oo ^ ' 



„2, 



here / (s) = is the density of W(l). 

The following is a well-known result 

P{|F(1)| > u\ 
lim — — = k 

for some k > 0; see, for example, Bertoin [9]. Hence, for fixed o < i, s 6 R, and as x — >• oo 

P{\Y(1)\ > x l l H (b - ay l l a \s\- l ' H } ~ - a)- 1 /"!^-!/^)- 

= kx- a / H {b-a)\s\ a l H . 
Now we apply the Dominated Convergence Theorem in equation (|4.30p to get 



\Z(b)-Z(a)\ > x} 



X— >oo 



lim L = k(b - a) / \s\°f"f B (s)d8 (4 32) 



= - a)2 a / 2 ^vr^ 1 / 2 r((a + H)/2H). 

Hence the constant in the theorem is C = k2 a l 2H -K^^T^a + H)/2H). □ 

Lemma 4.7. Let d = l ; 0<il<l and < a < 2. There exists a finite constant K > 
suc/i t/iat /or u > 1 , 

K -l u -a/H < p f gup > u \ < Ku~ a/H . (4.33) 

0<t<l J 
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Proof. Let S{t) = supo< s <j Then, by using the scaling property of W and condi- 

tioning, we have 



(4.34) 

Y ' 



sup \Z(t)\ > u } < F{ sup \W{x)\>u 
o<t<i J I \x\<S(l) 

= E(p{ sup \W(x)\ > 

It is well known that, for any e > 0, there exists a finite constant K such that for all u > 

p( sup |W(x)| > u\ < K exp ( - V (4.35) 
I \x\<i J V 2 + e/ 

See, for example, Lifshits [26j Section 14]. Consequently 

p{gP i | ZW |>,,}< K Eex P (- (2 + ^ (1)2 „ ). (4.36) 

Since, for all x > 0, the function ^(x) = exp ( — ^ 2 +^) x 2H ) ^ as P os itive derivative 

2#u 2 / u 2 \ 1 



2^V ( n 2 \ 

» W " 2 + e 6XP V (2 + e)x 2H J x^+i ' 



we derive 



sup \Z(t)\ > u 

0<t<l 

< Kn> jf" exp ( - ^P{S(1) > *} A, (4.37) 

where the last inequality follows from the change of variable x = u l / H y. Now by using the 
well-known estimate 

P{5(1) >y) < K(lAy~ a ), Vy>0, 
we obtain that for all u > 1, 

-a/H 



where .ftT > is a finite constant. This and (|4.37p together give the upper bound in (j4.33|) . 
The lower bound in (|4.33|) follows from Lemma 14.61 and the fact that 

P{|Z(1)| > u] = sup f{\Z(t)\ > u} < pj sup \Z(t)\ >u\. (4.38) 
0<t<l I 0<t<l J 
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The first equality in Equation (|4.38[) follows from the fact that the function 

t -> P{|Z(t)| >u}= P{|Z(1)| > t~ H/a u} 

is an increasing function for t £ (0, 1]. □ 

The following theorems are for laws of the iterated logarithm for the maximum local 
time L* ([r, r + h]) = sup xgR d L(x, [t,t + h]) and uniform Holder conditions of local times 
of a-time fractional Brownian motions. 

Theorem 4.2. Let d < a/H and suppose Y is not a subordinator. 

(1) . There exists a finite constant K > such that for any r > with probability 1 

L(x,t + h) — L(x,t) 

hmsup sup — —. - — — , : -, ; < K. (4.39) 

h _> F xe ^/ i 1 -^/«(loglog/i- 1 )^( 1 + 1 /«) 

(2) . For any T > 0, there exists a positive constant K such that almost surely 

L(x,t + h) — L(x,t) 
hmsup sup sup — : — ,„, , ; s < K. (4.40) 

Proof. Eq. (|4.39p follows from Lemma 14.71 and a chaining argument as that in the proof 
of Theorem 2 in [32] . The proof of Eq. (|4.40p , using Lemma 14.41 is very similar to that of 
Xiao |42} Theorem 3] and Ehm |19^ Theorem 2.1]. We omit the details. □ 

Theorem 4.3. Let d < a/H and suppose Y is a stable subordinator of index a < 1. 

(1) . There exists a finite constant K > such that for any r > with probability 1 

lim sup sup ^r + ^-^x r) 

(2) . For any T > 0, there exists a finite constant K > such that almost surely 

L(x,t + h) - L(x,t) 

lim SUp SUp SUp : ; r— 777/ < K- (4.42) 

The Holder conditions for the local time of a stochastic process X(t) are closely related 
to the irregularity of the sample paths of X(t) (cf. Berman [7]). In the following, we will 
apply Theorems 14.21 and 14.31 to derive results about the degree of oscillation of the sample 
paths of X(t). 

Theorem 4.4. Suppose Y is not a stable subordinator. Let X = {X(t), t € be an 

a-time fractional Brownian motion in M. d with FL < a. For any r £ R+, there exists a finite 
constant K > such that 

liminf sup „. ttt^t — ~ttwm i i ; ^ > K a.s. (4.43) 

^0 , e B(r,r-)^/«/(loglogl/r)^(l+V«) " V > 
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For any interval T C 



liminfinf sup „, — , X ^}\ : ; ; > K a.s. (4.44) 



/n particular, X(t) is almost surely nowhere differentiable in 



Proof. Clearly, it is sufficient to consider the case of d = 1, where the condition of Theorem 
g2] (i.e. 1 < a/H) is fulfilled. For any interval Q C R+, 

Ai(Q)= / L(x,g)dx<L*(g)f sup \X{s)-X{t)\ j. (4.45) 

Let Q = B(r,r). Then (fP3"|) follows immediately from ([439]) and ([OS]) . Similarly (fOI]) 
follows from (ET40]) and (H451 . □ 

Remark 4.2. Theorem\4-4\ extends partially the results obtained by Nane 130^ . 



Theorem 4.5. Suppose Y is a stable subordinator of index a < 1. Lei X = S R+} 

6e a-time fractional Brownian motion in M. d with H < a. For any r E R + , i/iere exists a 
finite constant K > swc/i i/iai 

liminf sup jy , >if a.s. (4.46) 

r^o seB (T,r)r H / a /(loglogl/r) H / a ~ 

For any interval T C R+ 

liminfinf sup - > if a.s. (4.47) 

r^o tzT sm l r) r H / a /(logl/r) H / a ~ 

In particular, X(t) is almost surely nowhere differentiable in R_|_. 

Acknowledgments. Authors would like to thank the two referees for their comments and 
corrections that helped improve the paper. 

5 Appendix 

As an appendix, we provide the following lemmas, which are used in the proofs of our main 
results in Section 2J Lemma 15.11 is from Xiao [42] , which is used to prove Lemma 15.21 



Lemma 5.1. Let < 7 < 1 be a constant. Then for any integer n > 1 and any x\, ■ ■ ■ ,x n £ 
R, we have 



1 1 



o min{\x -Xj\f,j = 1, • • • ,n} 
where K > is a finite constant depending only on 7. 



dx < K n 1 
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Lemma 5.2. Let < 7 < 1 be a constant. Then for any integer n > 1 and any x±, ■ ■ • , x n S 

Pi(x) 



J 

J e. 



dx<Kn^, (5.1) 



min{|x — Xj\i,j = 1, • • • , n} 
where K > is a finite constant depending only on 7 and a. 

Proof. We recall the following asymptotic bounds from |37] for the stable density function 
Pi(x) as x — > 00 (the asymptotics for the case x — > —00 are obtained by changing x to — x). 
For < a < 1: 

Pi 0*0 < ^^' 1+a ', as x — » 00. 
For a = 1 and /3 = (this is the symmetric Cauchy case): 

Pi(x) < K x~ 2 as x — > 00. 

For a > 1 and — 1 < f3 < 1: 

Pi(x) < isTx^ 1 " 1-0 ), as x — )• 00. 

For a > 1 and /3 = —1, 1: 

pi(x) < K max{x~( 1+Q ), x- 1+Q/2(Q - 1) exp(-c(a)x Q/(a - 1) )}, as x -> 00. 
Now we observe that the left-hand side of (15.111 can be written as 



W' +1 _ ?M , x 

— ^ / . mini t — t - 7 i = 1 ... r) V 



l min{|x - Xj|t, j = 1, • • • ,n} 

/M 
— : — n i — : ; r dx (5.2) 
-M nxm{\x-Xj\i,j = !,■■■ ,n} 

+ > max Di(x) / 



1 



|J|>Af 



min{|x + I — Xjpjj = 1, • • • , n} 



It can be verified that Equation (|5.ip follows from (|5.2p . the asymptotics of pi(x) and 
Lemma 15.11 □ 



Lemma 15.31 is taken from Ehm [19] and Lemma 15.41 is due to Cuzick and DuPreez [14 
(the current form is from Khoshnevisan and Xiao 



Lemma 5.3. For any integer n > 1, and /3j £ (0, 1) for 1 < j < n, for all h > 0, we /iai>e 

f n 1 

= h n-^ n^rq-/^) 

r(i + n-E?=i^)' 



24 



Lemma 5.4. Let £1, • • • , £ n be mean zero Gaussian variables which are linearly independent, 
then for any nonnegative function g : R — > R+, 



g(v!)exp 



1 / - 
-Var J2 V & 



j'=i 
(2vr)(™-i)/2 



(detCov(£i 



d«i • • • dv n 



00 / v 

-oo V 0"l 



(5.3) 



where detCov(£i, • • • , £ n ) denotes the determinant of the covariance matrix of the Gaussian 
random vector (£1, . . . , £ n ), and where a\ = Var(£i|£2, ■ ■ ■ , £n) *s i/ie conditional variance of 
£1 ^roen £ 2 , • • • ,£n- 
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